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Abstract 

A single gene, regulating its own expression via a positive feedback loop, constitutes 
a common motif in gene regulatory networks and signalling cascades. Recent experi- 
ments on the development of competence in the bacterial population B. suhtilis show 
that the autoregulatory genetic module by itself can give rise to two types of cellular 
states. The states correspond to the low and high expression states of the master 
regulator ComK. The high expression state is attained when the ComK protein level 
exceeds a threshold value leading to a full activation of the autostimulatory loop. 
Stochasticity in gene expression drives the transitions between the two stable states. 
In this paper, we explain the appearance of bimodal protein distributions in B. suh- 
tilis cell population in the framework of three possible scenarios. In two of the cases, 
bistability provides the basis for binary gene expression. In the third case, the system 
is monostable in a deterministic description and stochasticity in gene expression is 
solely responsible for the appearance of the two expression states. 

1 Introduction 

Positive feedback loops are common motifs in gene transcription regulatory networks and 
signaling cascades. The simplest such motif is the autoregulatory loop in which the pro- 
teins synthesized by a gene stimulate the production of more proteins in an autocatalytic 
fashion [HE], [3]. In most cases, the concerned gene is also expressed at a basal level, i.e., 
proteins are synthesized even when the positive feedback is non-functional. The autoregu- 
latory dynamics have a nonlinear character and this combined with positive feedback may 
give rise to binary gene expression in a range of parameter values. The protein levels, as a 
result, have a bimodal distribution in a population of cells. In a fraction of cells, the protein 
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level is low and in the rest of the population the level is high. Theoretical predictions of 
binary gene expression have been verified in experiments on single gene autocatalytic mod- 
ules in bacteriophage A and S. cerevisiae ^ |4j . The experimental findings further suggest 
that the observed bimodality has a stochastic origin. A remarkable example of population 
heterogeneity, brought about by a combination of autoregulatory positive feedback and 
stochasticity, is provided by the bacterial population B. subtilis in which a fraction of the 
population develops genetic competence. Microorganisms like bacteria have to cope with a 
multitude of antagonistic agents and environmental conditions in order to live. Under such 
circumstances, the bacteria may adopt a number of strategies to optimize their chances 
of survival [H [5]. One such strategy is the development of genetic competence, observed 
in some bacterial organisms. In the competence state, specialized proteins are synthesized 
which allow the cell to take up large pieces of DNA from the environment and incorporate 
them into the bacterial genome. New traits are thus acquired from genetically distinct or- 
ganisms. Experiments show that only a small fraction of the bacterial population reaches 
the competence state. The resulting phenotypic diversity in the population may prove to 
be advantageous. The individual cells in a homogeneous population share the same fate 
when subjected to harmful infiuences. Diversity enhances the chance that a fraction of the 
population, even if small, is able to survive and adapt to the changed circumstances. In 
B. subtilis, the development of competence is regulated by the transcription factor ComK 
synthesized by the comK gene. The protein functions as a master regulator which acti- 
vates the transcription of several genes including those necessary for DNA uptake. The 
ComK activity in turn is controlled by a host of other proteins. An autoregulatory posi- 
tive feedback module forms the core of the complex regulatory network. ComK binds to 
the promoter of its gene and promotes its own production. The positive feedback gives 
rise to bimodality in the cell population with low and high comK expression states as 
the stable states. In the competence state, the level of ComK proteins is high enabling 
ComK to act as a transcription factor. Two independent experiments [6l [7] have confirmed 
that an autostimulatory loop of comK expression is by itself sufficient to establish compe- 
tence bimodality in a bacterial culture. The experimental findings moreover suggest that 
stochasticity plays an essential role in the establishment of competence. 

In this paper, we study a simple model of autoregulatory positive feedback involving 
a single gene using both deterministic and stochastic descriptions. In the deterministic 
case, positive feedback and nonlinearity result in bistability in a range of parameter val- 
ues. The two stable steady states correspond to low and high gene expression levels. The 
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high expression state is reached when the protein level exceeds a threshold value. The 
autocatalytic switch is then triggered bringing about a full activation of the autostimula- 
tory loop. In the absence of such activation the proteins are synthesized at a low level. 
Bimodality in a cell population requires the autocatalytic switch to be triggered in a frac- 
tion of the cell population. This is where stochasticity in gene expression comes into the 
picture. Several recent studies, both theoretical and experimental, highlight the significant 
role of stochasticity in gene expression and its regulation [HI [9l [TOl [HI, [12] . The two stable 
steady states are separated by an unstable steady state. The corresponding protein level 
(intermediate level of gene expression) provides the threshold for the triggering of the au- 
tocatalytic switch. The low (high) expression state is obtained when the protein level is 
below (above) the threshold value. Stochasticity in gene expression gives rise to fluctua- 
tions in the protein levels and the fluctuations, if sufficiently large, bring about transitions 
across the threshold. In the deterministic picture, bifurcations occur at two special values 
of the parameter Jq, the rate for basal protein synthesis. At the lower (upper) bifurcation 
point, there is a transition from monostability (bistability) to bistability (monostability). 
This framework provides an alternative explanation of population heterogenity. Inducer 
molecules are often required to initiate gene expression at the basal level. The distribution 
of the molecules may be non-uniform in a population of cells. Thus, the basal levels in 
the individual cells are not identical but have a disribution around an average value. If 
this distribution overlaps with the upper bifurcation point, the cell population develops a 
bimodal character. There is also a third explanation for population heterogeneity which 
is solely based on stochasticity in gene expression. In this case, the system is not bistable 
in the deterministic picture and bimodality occurs due to random transitions between the 
low and high expression states. In this paper, we explore the basis of bimodal protein 
distributions in the three scenarios outlined above. The results are interpreted in terms of 
the development of genetic competence in B. subtilis bacterial population. 

2 Deterministic model 

We consider a simple model of autoregulatory gene expression involving a single gene. The 
proteins synthesized by the gene form dimers. The dimer molecules bind to the promoter 
region of the gene and activate gene expression, thus constituting a positive feedback loop. 
Apart from autoactivation, the gene synthesizes proteins at a basal level. The detailed 
kinetic scheme of the model is shown in figure 1(a). The gene can be in two possible 



3 




FIG. 1(a): The kinetic scheme describing autoregulatory gene expression. G and G'*are 
the inactive and active states of tlie gene. In the inactive state, proteins are synthesized at 
a basal rate Jq. The protein molecules form dimers P2 which bind to the promoter region 
of the gene and activate the state G to G*. 




FIG. 1(b): The reduced kinetic scheme with effective activation and inactivation rate 
constants ^^(a;) and k^. 



4 



states G (inactive) and G* (active) . In the active state proteins are synthesized with rate 
constant Ji. In the inactive state G, "leaky" gene expression occurs at the basal rate Jq 
(Ji >> Jo)- The basal rate may be enhanced using appropriate inducer molecules. The 
synthesized proteins dimerize with K being the equilibrium dissociation constant. The 
protein dimer P2 binds to the gene in its inactive state G and activates the gene to the 
state G* . The rate constants ka and kd are the activation and deactivation rate constants. 
The synthesized proteins are degraded with a rate constant kp. The kinetic scheme in figure 
1(a) can be mapped onto a simpler scheme shown in figure 1(b). The eff'ective activation 
and deactivation rate constants k!^{x) and k!^ are given by 

1 ' / \ 1 1 ' 1 / \ 

^'^^'^^ = ^"l + (x/fc,)2' = 



where x denotes the protein concentration and ks = ^^K. 

In the simplified kinetic scheme of figure 1(b), the rate of change of protein concentration 
is given by 

dx ^ Jik'^jx) _^ Jpkd _^ ^ ,^ 
dt k'^{x) + kd k'^{x) + kd ^ 



In the steady state, = and one can identify a parameter region in which the system is 
bistable, i.e., has two stable steady states. These states correspond to low and high values 
of X. An unstable steady state (intermediate value of x) separates the two stable steady 
states. Figure 2 shows a plot of versus Jo where denotes the steady state protein 
concentration. The solid branches represent stable steady states and the dotted branch, 
the unstable steady states. In a range of parameter (Jq) values, the system is bistable. 
The other parameters have values ka = 0.0008, kd = 0.0005, kg = 500.0, Ji = 0.1 and 
kp = 0.0001 in appropriate units. Bistability is, in general, accompanied by hysteresis 
[HI [IB]. Let us assume that the system is in the lower steady state and the value of Jo is 
small. As Jq is increased (say, with the help of inducer molecules), the system continues to 
be in the low expression state. At a critical value Jquc, a discontinuous transition to the 
upper stable steady state occurs. If Jq is increased further, the system is monostable, i.e., 
there is only one stable steady state (the upper state). If the value of Jq is now reduced 
below Jquc, the system remains in the upper steady state which is a hallmark of hysteresis. 
At a lower critical value of Jq = Jolc (marked by a vertical line on the horizontal axis of 
figure 2), a transition from the upper to the lower stable steady state occurs. 

Hysteresis promotes robustness as once the system is in the upper stable steady state. 
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FIG. 2: Bistability and hysteresis, the solid (dotted) lines represent stable (unstable) steady 
states; x'^ is the steady state concentration of proteins and Jq, the basal rate of protein 
synthesis, serves as the bifurcation parameter. The short vertical lines on the horizontal 
axis denote the lower and upper bifurcation points, Jqlc and Jquc- 

small fluctuations of Jq around Jquc will not give rise to a transition to the lower stable 
steady state. Let us now assume that the basal gene expression is initiated with the 
help of inducer molecules. These molecules may have a heterogeneous distribution in the 
cell population (each individual cell contains the autoregulatory module) which gives rise 
to a distribution in the basal rates Jq. If the threshold value Jquc falls within the Jq 
distribution, the cell population exhibits bimodality. Cells in which the basal rate Jo is 
less (greater) than Jquc, are in the low (upper) stable steady state. Figures 3(a) and (b) 
illustrate this for a normal distribution of basal rates with mean = 0.00445 and variance 
= 0.0005. In figure 3(b), p{x) describes the steady state distribution in the protein levels. 
In the steady state, ^ = in equation ([2]), from which the basal protein synthesis rate Jo 
can be expressed as a function of x i.e., Jq = f{x). Let p(Jo) be the distribution in basal 
levels {Jo/kp is the steady state basal level). One can then write 

Pix) = p{jo) \jo=f(x) 1^1 (3) 

This way of explaining bimodality is consistent with an earlier proposal on the origin 
of binary gene expression [14]. We now discuss the other two mechanisms for obtaining 
bimodality taking stochasticity in gene expression explicitly into account. 
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FIG. 3(a): Normal distribution describing heterogeneous inducer distribution overlaps with 
the upper bifurcation point Jquc- 
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FIG. 3(b): The bimodal distribution in protein levels due to the heterogeneous distribution 
of inducer molecules. 
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3 Stochastic origins of bimodality 



We consider a simple stochastic model corresponding to the kinetic scheme in figure 1(b). 
In the model, the only stochasticity arises from the random transitions of the gene between 
the inactive and active states as in the minimal model of Cook et al. fl5|. Protein synthesis 
from the inactive (basal expression) and active states of the gene and protein degradation 
occur in a deterministic manner. We would like to determine the distribution of protein 
levels in the steady state of the cell population. Following the method outlined in [H], the 
concentration of proteins evolves as 

dx 

— = Jiz + Jq{1 - z) - kpX = f{x,z) (4) 

where z = 1(0) when the gene is in the active, G* (inactive, G) state. The random variable 
z switches values with stochastic rate constants k'^{x) (0 1) and k'^ (1 0). Let pj{x,t) 
(j = 0, 1) be the probability density function when z = j. The total probability density 
function is 

p{x,t) = po{x,t) + pi{x,t) (5) 
The rate of change of probability density is given by 

= -§^[f{x,j)PjM]+^[W,,p,{x,t) - W,,p,{x,t)] (6) 

where Wkj is the transition rate from the state k to the state j and Wjk is the same for the 
reverse transition. The first term in equation ([6]) is the so called "transport" term repre- 
senting the net fiow of the probability density. The second term represents the gain/loss in 
the probability density due to random transitions between the state j and other accessible 
states. In the present case, equation (I6|) gives rise to the following two equations: 

= ~^{(^o - kpx)p{x,t)]} + kdPi{x,t) - k'^{x) Po{x,t) (7) 

= "^{(^1 ~ ^p^) + k'ai^) Po{x,t) - kdPi{x,t) (8) 

Using equation ([5]) , the steady state solution of equations (jTj) and (jHD is given by 

p{x) = C{kpX- Jo^iJi - kp x)~^+^(x2 + k^J'^Explu] (9) 
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where 

ka arctan(a;//Cs) Jq ks Jl{—ka + kp) ~\~ kp kg /lq, kp kg , , 

Jq + fcp /c^ ^p(<^o + ^s) 2(Jo + fcp fc^) 

and C is the normalization constant. 

Figures 4(a)-(d) show the plots of versus x as the activation and deactivation rate 
constants ka = 0.0008*/;, and k^ = 0.0005*/;. are progressively changed by varying the factor 
h. The other parameters have values kg = 500, Jq = 0.0035, Ji = 0.1 and kp = 0.0001. As h 
is changed from 2 to 100, there is a transition from unimodality (figure 4(a)) to bimodality 
(figures 4(b) and (c)) to again unimodality (figure 4(d)). The unimodal distributions 
correspond to low (figure 4(a)) and high (figure 4(d)) gene expression states. In all the 
four cases, the deterministic dynamics (equation ((21)) lead to bistability in the steady state 
(figure 2). The three steady state solutions for Jq = 0.0035 are XstaUei = 50.11, XunstaUe = 
160.60 and Xstabie2 = 418.13. The bimodality observed in figures 4(b) and (c) are due to 
stochastic transitions between the stable steady states brought about by the fiuctuations 
associated with the protein levels. The two stable steady states XstaUei and Xstabie2 are 
separated by the unstable steady state Xunstabie- When protein levels are below (above) 
^unstable, the low (high) cxprcssiou state becomes the stable steady state. Fluctuations in 
the protein levels (due to stochastic gene expression) are responsible for the excursions 
from one state to the other. In the deterministic picture (figure 2), the system is in the 
lower stable steady state for the parameter value Jq = 0.0035 used in obtaining the plots 
in figure 4. Figures 4(b)-(d) thus clearly demonstrate that noise can alter the deterministic 
outcome in a significant manner. In case (d), the fiuctuations associated with the lower 
protein level are so strong that excursions to the higher protein level occur with probability 
one. In terms of the autoregulatory genetic module, the autostimulatory feedback loop is 
fully activated when the protein level x is > Xunstabie so that the high expression state is 
achieved in the steady state. In the deterministic picture, the time evolution of a dynamical 
system can be predicted with absolute certainty once the parameter values and the initial 
state are specified. In the present case, the different rate and binding constants constitute 
the parameters. The state of the system at time t is given by the amount of proteins x{t). 
The value of x{t) is obtained by solving the differential equation (equation (2)) for a fixed 
set of parameter values and with a knowledge of the initial state x{to) at time to. The 
time evolution of the system is represented by a trajectory in state space (one-dimensional 
in the present case). The trajectory starts from the point x(to)and ends at a fixed point 
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= 0) describing a stable steady state. In the region of bistability, the two stable steady 
states Xstabiei and Xstabie2 have their individual basins of attraction [l7l[IH]- A trajectory 
which starts in one particular basin of attraction reaches the corresponding stable steady 
state in the course of time. The time evolution of a system stops once the steady state is 
reached. A steady state is stable (unstable) if the system comes back to it after a weak 
perturbation is applied. Small fluctuations in the protein level x leave the system in the 
same basin of attraction. There may, however, be excursions from one basin to the other 
when the fluctuations are of sufficiently large magnitude. The probability of transition 
from one basin of attraction to the other depends amongst other factors on the value of 
Jo, the basal rate of protein synthesis. The gap between XunstaUe and XstaUei is smaller and 
that between Xstabie2 and XunstaUe larger as Jq approaches Jqc/c- The reverse situation is 
true as Jq approaches the lower bifurcation point. The plots in figure 4 have been obtained 
for progressively higher values of the activation rate constant ka- The value of Jq = 0.0035 
is closer to the upper bifurcation point Jquc- The protein fiuctuations are amplified for 
higher values of ka- The fiuctuations have to bridge a smaller gap for transition from the 
basin of attraction of XstaUei'^o that of Xstabie2 than in the case of the reverse transition. In 
the case of figure 4(a), the system remains in the basin of attraction of XstaUei- As ka is 
made higher, a greater fraction of the cell population attains the high expression state. In 
the case of figure 4(d), almost the whole cell population is in the high expression state. 
Figures 5(a)-(d) show plots similar to those in figure 4 for a lower value of Jq = 0.0030. 
The gap between Xunstabie and Xstabiei is now larger and that between Xstabie2 and XunstaUe 
smaller (see figure 2). The balance in this case tilts in the favour of the lower stable steady 
state. 

We now consider the third case in which a bimodal protein distribution has a purely 
stochastic origin. The system is monostable in the deterministic description. An earlier 
study by Kepler and Elston [Hj provides examples of such cases. Some other studies 
have explored the basis of stochastic binary gene expression in different settings (without 
positive feedback) [HI [121 [201 El] ■ An example in the case of autoregulated gene expression is 
shown in figure 6 for the parameter values ka = 0.0012, k^ = 0.0004, ks = 500.0, Jq = 0.01, 
Ji = 0.1 and kp = 0.0001. In the deterministic description, there is only one stable steady 
state, x^ = 581.3. The protein distribution is obtained from the analytic expression given 
in equation (9). The stochastic model considered in this section is analytically tractable 
because of certain simple assumptions. The only stochasticity considered in the model is 
that associated with random gene activation and deactivation. The autocatalytic feedback 
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FIG. 4: The steady state distribution, p(a;)versus x, in protein levels with the activation 
and inactivation rate constants given by ka = 0.0008 */i and ka = 0.0005 The plots are 
obtained for different values of h, (a) h = 2, (h) h = 20, (c) h = 50 and (d) h = 100. The 
basal rate of protein synthesis is Jq = 0.0035. The three steady states in the deterministic 
case are Xgtabiei = 50.11, x^nstaUe = 160.6 and Xstahie2 = 418.13. 
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FIG. 5: The steady state distribution, p{x) versus x, in protein levels. The parameter Jq 
has the value Jq = 0.0030 with ka — 0.0008 * h and kd — 0.0005 * /i as in the case of figure 
4. The values of h are (a) /i = 2, (b) /i = 20, (c) h ^ 50 and (d) h = 250. The three steady 
states in the deterministic case are XstaUei — 39.55, XunstaUe — 178.22 and Xstabie2 — 409.14. 
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is incorporated in an effective rate constant k'^{x). The gene expression is considered 
as a one-step process, i.e., the intermediate stage of mRNA synthesis is not explicitly 
taken into account. We now describe the results of a detailed simulation based on the 
Gillespie algorithm [22] which takes the two-step nature of gene expression into account and 
treats the distinct biochemical events to be stochastic in nature. The different biochemical 
reactions are listed in equations f|TTll-(l2Tll: 



G + P2^GP2 (11) 

GP2 -^G + P2 (12) 

GP2^G* (13) 

G* GP2 (14) 

G*^m (15) 

G^m (16) 

m — i> deg. (17) 

m^P (18) 

P deg. (19) 

P + P^ P2 (20) 

P2^P + P (21) 

In the above equations, the mRNA and protein are represented by m and P respectively, P2 
is a protein dimer and GP2 denotes the intermediate state of a protein dimer bound to the 
gene in its inactive state G. Equations (fTTll and ( fT9l ) describe the degradation of the mRNA 
and protein molecules. In the simulation, the stochastic rate constants associated with the 
equations ([II|)-(l2ll) are c(l) = 0.003, c(2) = 0.16, c(3) = 0.004 * h, c(4) = 0.0006 * h, 
c(5) = 0.1, c(6) = 0.0015, c(7) = 0.0001, c(8) = 0.000001, c(9) = 0.008, c(10) = 0.01 and 
c(ll) = 0.01 in appropriate units. The simulation is carried out for three different values 
of h = 1, 10 and 200 respectively. The results are shown in figures 7(a)- (c). The plots 
on the left show the time trajectories, x(t) versus t, where x(t) is the amount of proteins 
at time t. The plots on the right show the distributions p{x) versus x on repeating the 
simulation 3000 times. The quantity p{x)dx provides a measure of the fraction of cells in a 
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FIG. 6: The steady state distribution, p{x) versus x, in protein levels for ka = 0.0008, 
kd = 0.0004 and Jq = 0.001. The bimodal distribution has a purely stochastic origin. In 
the deterministic case, there is only one stable steady state = 581.3. 

population with protein levels between x and x + dx. For the parameter values used in the 
simulation, only a small fraction of cells is in the high expression state. Figures 7(a) and 
7(b) show that the stochastic nature of the biochemical events involved in gene expression 
is responsible for a bimodal protein distribution. 

4 Discussions 

In this paper, we study how positive feedback combined with stochasticity gives rise to 
binary gene expression, i.e., a bimodal distribution in the protein levels in a population of 
cells. There are some earlier studies [21 [19] on the same issue but the modeling detail and 
context are different. The motivation for the present study comes from the experimental 
observation that a single comK gene, which autoregulates its expresssion via a positive 
feedback loop, is by itself sufficient to generate heterogeneity in a population of B. suhtilis 
[SI El [7]. A fraction of the cell population develops competence due to the high expression 
state of comK. This is so when the ComK protein level exceeds a threshold value thus trig- 
gering the full activation of the autostimulatory loop. The comK autoregulatory genetic 
module is at the core of a complex network of molecular interactions which regulate comK 
transcription and the stability of the ComK proteins. In this case, only a small fraction 
of cells, about ten percent, develops competence. When experiments are carried out on 
the isolated genetic module, the fraction of cell population in the high comK expression 
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FIG. 7: Results of simulation based on the Gillespie algorithm. The stochastic rate con- 
stants for gene activation and inactivation are c(3) = 0.004 * h and c(4) = 0.0006 * h The 
factor h has values (a) /i = 1, (b) /i = 10 and (c) h = 200. The plots in the first column 
show the variation of protein amount as a function of time. The plots in the second column 
show the corresponding steady state distributions, p{x) versus x. 
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state can be quite large. It has further been suggested that stochasticity in gene expres- 
sion is responsible for throwing the autocatalytic switch [S, El [7]. In our study, we focus 
only on the core module, namely, the comK autoregulatory genetic module of the net- 
work regulating competence development. The module represents a single gene (comK), 
the protein product of which autoactivates its own synthesis through dimerization and 
subsequent binding at the appropriate region of the DNA. Our model incorporates these 
minimalist features of the autoregulatory comK module. We have explored the basis of 
binary gene expression in the framework of three different possibilities. In the first two 
cases, the underlying dynamics lead to bistability in a deterministic description. A hetero- 
geneous distribution in inducer molecules may give rise to a bimodal distribution in the 
protein levels. In the second and third cases, we take stochasticity in gene expression into 
account and derive an analytic expression (equation Q) for the steady state distribution of 
protein levels. The analytical tractability of the stochastic model arises from two assump- 
tions. Firstly, the two major steps of gene expression, namely, transcription (synthesis of 
mRNAs) and translation (synthesis of proteins) are combined into a single step leading 
to protein production. Secondly, the only source of stochasticity in the model lies in the 
random activation and deactivation of the target gene expression. The first assumption 
provides the basis for several studies of stochastic gene expression [I5l [HI ED, [23]. The 
second assumption is strictly valid when the dominant source of noise is associated with the 
random activation and deactivation of gene expression. This is so in the case of slow pro- 
moter kinetics. As discussed in detail in [12], slow transitions between the promoter states 
result in transcriptional bursts of mRNA synthesis and increased heterogeneity within a 
cell population including bimodal protein distributions. Experimental evidence of tran- 
scriptional bursting has been obtained for both prokaryotes and eukaryotes [241 ESj. A 
recent experiment on stochastic mRNA synthesis in mammalian cells [26] shows that the 
mRNA levels display large cell-to-cell variations due to random, infrequent activation of 
gene expression. The statistics of the variations are adequately described by a model in 
which the only source of stochasticity lies in the random activation and deactivation of the 
gene. There could be a number of factors which lead to slow transitions between the pro- 
moter states. Chromatin remodeling has been conjectured to cause transcriptional bursts 
in eukaryotic systems [12]. In both prokaryotes and eukaryotes, pulsatile gene expression 
may also result from regulatory molecules binding at and unbinding from the DNA sites, 
DNA undergoing conformational changes so that the RNA polymerase has only brief ac- 
cess to the promoter region etc. In the case of E. coli, there is experimental evidences 
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that long periods of inactivity are interspersed by shorter periods when the gene is in the 
active state [24]. An earlier result of Ozbudak et al. on B. subtilis [8] has been reanalyzed 
to show that the data are not inconsistent with the possibility of transcriptional bursts 
(random gene activation-inactivation) (24]. Our simple stochastic model, based on random 
transitions between inactive and active gene states, is thus consistent with experimental 
reality. In the case of bistability, stochasticity triggers transitions between the two stable 
steady states which is responsible for bimodal protein distributions. In the case when the 
system is monostable in a deterministic description, binary gene expression can still occur 
due to a combination of positive feedback and stochastic transitions between the inactive 
and active states of the gene. Our analytical results are supported by the simulation results 
in the case of a more detailed stochastic model in which transcription and translation are 
treated as separate processes and stochasticity associated with all the biochemical steps 
(equations (fTTl)-(l2Tl)) are taken into account. In both the cases, the results are valid over 
a wide range of parameter values. We now briefly discuss the experimental possibility for 
distinguishing between the three mechanisms discussed in the paper. As shown in figure 
2, bistability implies hysteresis. A properly designed experiment can detect hysteresis in 
the response {x'^ in figure 2) as the variable along the x — axis (Jq, the basal rate of pro- 
tein synthesis in figure 2) is changed. Discontinuous jumps in response at the bifurcation 
points and a non-reversible response are the hallmarks of hysteresis. The value of Jo may 
be changed using appropriate inducer molecules. One can use a cell sorter and separate a 
subpopulation from a bimodal cell population. The subpopulation develops bimodality in 
the course of time if there are stochastic transitions between the low and high expression 
levels. 

Siiel et al. [2^ have investigated competence development on the basis of a model 
describing an excitable stochastic system. The key ingredients of the model are: the 
comK autoregulatory loop, the inhibition of ComK degradation by ComS proteins and 
repression of the comS gene by ComK. Theoretical analysis of the model dynamics is 
combined with experiments to gain insight on the entry into and exit from the competence 
state. This state corresponds to an unstable fixed point of the model dynamics. The 
system has only one stable steady state in which the ComK level is low. Fluctuations 
in the levels of ComK/ComS excite the system into the competence state with eventual 
return to the noncompetence state. In the excitable system, repeated stochastic triggering 
of the competence state is thus possible. Some of the premises of the model like the 
"indirect" repression of the comS gene by ComK need experimental confirmation under 
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wild-type expression conditions [28]. The study nonetheless is an elegant example of how 
model studies combined with experiments can provide a new perspective on noise-induced 
phenomena in biological systems. Our model has the comK autoregulatory loop as the 
sole ingredient and focuses on the specific experiment by Smits et al. p] on the single 
autoregulatory module. Genetic competence in B. subtilis provides a concrete example of 
a natural system in which a single gene, regulating its expression via an autoregulatory 
positive feedback loop, is by itself sufficient to establish two types of stable states in the cell 
population. Recently, two groups have independently discovered a similar phenomenon in 
the human fungal pathogen Candida albicans [2HI30]. In both the cases, the autoregulatory 
modules are parts of complex genetic circuitry. The single gene modules almost exclusively 
control the cellular switch operating between two stable states. The resulting heterogeneity 
is epigenetic in nature. B. subtilis and Candida albicans thus illustrate the essentiality and 
sufficiency of network modules in explaining particular types of biological function. The 
role of the other components of the associated regulatory networks lies in modulating the 
functional response. In B. subtilis, several genes regulate the expression of the comK gene 
the protein product of which regulates the expression of several other genes. The products 
of the regulatory genes modulate the threshold for the triggering of the autocatalytic 
switch and influence the stability of the ComK proteins. The additional circuitry probably 
includes features which further stabilize the steady states. In Candida albicans, the WORl 
gene acts as the master regulator. The gene autoregulates its own expression via a positive 
feedback loop. The switch now operates between the cellular states: white and opaque. 
The two types of cells, white and opaque, differ in their morphologies, the genes they 
express, the host tissues in which they are resident and also in their mating characteristics. 
In the white cells, WORl is expressed at low levels whereas the levels are high in the 
opaque state. As in the case of B. subtilis, stochasticity appears to drive the transitions 
between the two types of cell. The results derived in this paper, specially those pertaining 
to the combined effects of bistability and stochasticity, should be of relevance in explaining 
the white-opaque switching in Candida albicans. 
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